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In the present paper, we develop in detail a fully geometrical method for deriving perturbation 
equations about a spatially homogeneous background. This method relies on the 3 + 1 splitting of 
the background space-time and does not use any particular set of coordinates: it is implemented in 
terms of geometrical quantities only, using the tensor algebra package xTensor in the xAct distri- 
bution along with the extension for perturbations xPert. Our algorithm allows one to obtain the 
perturbation equations for all types of homogeneous cosmologies, up to any order and in all possible 
gauges. As applications, we recover the well-known perturbed Einstein equations for Friedmann- 
Lemaitre-Robertson-Walker cosmologies up to second order and for Bianchi I cosmologies at first 
I order. This work opens the door to the study of these models at higher order and to that of any 

^ . other perturbed Bianchi cosmologies, by circumventing the usually too cumbersome derivation of 

the perturbed equations. 
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INTRODUCTION 



o 

u 

^ . . . 

• Cosmological perturbation theory constitutes the cornerstone of our current understanding of the origin, evolution 
I I and formation of large-scale structures. The evolution history of perturbations is written about a fixed homogeneous 
O . and isotropic Friedmann-Lemaitre-Robertson- Walker (FLRW) background space-time, and the interpretation of cos- 
> ' mological observations (such as WMAP, QUIET and ACT within this model, converge toward a coherent and 

^ . unified picture of the Universe. This picture is likely to get even clearer when the next generation of large-array cosmo- 
I— I' logical observations (such as EUCLID and SKA), becomes operational. It is widely expected that these observations 

will generate large amount of data that will provide a percent level accuracy for the cosmological parameters. 
^ ' Cosmological perturbation equations are simple and straightforward to derive at linear order, but they are inade- 
, quale for understanding the late-time evolution of the Universe, precisely when the non-linear gravitational effects, 
' carrying the information of the physics of current interest, come into play. Going beyond first order is a difficult task, 
^-H , and in some cases it becomes extremely arduous to even perform a coordinate or gauge transformation at non-linear 
^^O • order by hand. To the best of our knowledge, there is no available easy-to-use software designed for cosmology and 
I capable of deriving all equations of motion for perturbed variables. The closely related available option is the GRTen- 
■ sor package Q, which runs on Maple or Mathematica. However, the outputs generated at linear order are already a 
' bit complicated to understand, let alone its outputs at non-linear orders, as it relies exclusively on a properly defined 
set of background coordinates each time it acts on a perturbed variable. 

To fill up this gap, we have developed an algebra package for cosmological perturbation theory, called xPand [Bj, 
which uses the tools of the tensor algebra package xTensor and an extension for perturbation, xPert The xTensor 

^ . and xPert are part of the xAct distribution [g that run on Mathematica and they are available under the General 
d , Public License. xPert is specifically designed to perform perturbations on arbitrary background space-times 0, Q, 
but it lacks the features for specializing to a specific background space-time and a specific form for the metric 
perturbations, as is needed in the case of cosmology. In [^,[1^ these packages were used to study perturbations about 
a spherically symmetric space-time, more precisely around a Schwarzschild solution of the Einstein field equations. 
The xPand package now allows one to derive, in a simple and user-friendly manner, all the necessary equations for 
cosmological perturbation variables, around any homogeneous background space-times, at any order and in any gauge. 
Specifically, the available type of cosmologies cover the Minkowski, FLRW (fiat and curved), and Bianchi space-times, 
and the available predefined gauge are: general gauge (no gauge choice), comoving gauge, flat gauge, isodensity gauge, 
Newtonian gauge and synchronous gauge. 

The paper is organized as follows. In section HI we provide a general overview of the mathematical framework on 
which xPertis built. In sectionllll we detail the 3 + 1 splitting of the background manifold into a family of homogeneous 
hypersurfaces orthogonal to a fundamental observer's velocity. In section lllll we decompose the perturbed metric with 
respect to this foliation and define the scalar, vector and tensor perturbations. Each of these sections is supplemented 
by a presentation of the associated implementation in xPand, by means of several detailed examples. Finally, we 
summarize and discuss in section HVl the features and performances of the package. 
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I. PERTURBATIONS AROUND A GENERAL SPACE-TIME 

In this section we briefly review the algorithm of xPert Q, which constitutes the basis of our method. For more 
details about perturbation theory in the context of cosmology, we refer the reader to, e.g., [13, E^- 



A. General framework 



Let us consider in what follows a background manifold Ai along with its perturbed manifold Ai . Both are related 
by means of a diffeomorphism (j) : M — > A^. Tensorial quantities are thus transported from one manifold to the other 
with the help of the pull-back (/)*, the push- forward 0^,, and their respective inverses. The metric of the perturbed 
manifold relates to that of the background as 

°° A"fnl 

0'^(g)=g + A[g]=g + ^ Ji. (1) 

Tl=l 

Here and in the sequel, we use boldface symbols for tensorial quantities, an over-bar for background quantitie^, and 
the notation A[T'] (resp. A"[T']) for the total (resp. n"^ order) perturbation of a tensor T. One may prefer to write 
instead of A, as the definition of the perturbations depends on the diffeomorphism (p, that is on the choice of 
the gauge. We however choose to omit this reference for the sake of clarity, and in order not to burden the notation 
unnecessarily, we shall moreover use the short-hand: T = (f)*{T), for any perturbed quantity. 

Unless otherwise specified, when we write down the components of a tensor, these should be understood as expressed 
in a general basi^ (this holds equally for the background and perturbed tensors, and for the perturbations). Since 
all perturbation orders live on the background manifold, as they are the result of the pull-back of a tensorial quantity 
living on the perturbed manifold, we shall raise and lower indices using the background metric. We have for instance: 

(">;i'^-=5A-^5-("V, (2) 
for the n*^ order metric perturbations ^"^h = A"[g]. 



B. Expansion of the curvature tensors 



1. Mathematical framework 

The inverse of the metric tensor is given by the relation: 

g-' = {g + A[g])-' . (3) 

Expanding it into 

oo 

and making use of the definition ^"^ft, = A"[g], we obtain the n"^ order perturbation of g^ : 

i\ 

fci! . . . km^. 



A" 



where the sum X](A: ) I'^ns over the 2" ^ sorted partitions of n for m < n positive integers, such that fci + . . . km ~ n. 
Note that A" [(.g~^)'"'] ^ = .g^'^.g'"" (c-g-, we have at first order: I^'^[{g-'^f\ = - ^^^hf""). 



^ This convention differs from the one adopted in Q. We however opt for this choice as it reflects the standard usage in cosmological 

perturbation theory. 
^ We use Greek letters (o, /3, /x, v, p, ...) for space-time indices. 
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By means of relation ([5]) , we can express the perturbation of the connection components as Q 

where the last term of the right-hand side is defined by 

= \ (V, <">/i^^ + ^" V ' ^">/v.) . (7) 
The perturbation of the Riemann tensor is given in ah generahty by 

Tl-l ^ \ 

(8) 



(9) 









"P. 



k=l 

and for a connection compatible with the metric, we have, from equation 

m -v 



The symbol (/.i o t^) denotes the repetition of the preceding expression with indices /i and v exchanged. The 
perturbation of the Ricci tensor is simply obtained by contracting the second and fourth indices of A" [-R^i/p*^] in the 
previous expression, and the perturbation of the Ricci scalar, R = gP°' Rp^, is written: 

At last, the perturbation of the Einstein tensor is expressed according to 

^ n k I 

A" [G,.] = A" [i?,.] - 2 E E fc,,,fa"\-_fcv V [5^1 A"-^-'^- [i?,.] . (11) 

2. Implementation in xPert 

All the perturbative expansions expounded above are implemented in the package xPert 0- For completeness, we 
here briefly review its main commands. The package can be loaded by evaluating: 

In[l] := «xAct'xPert' 

(Version and copyright messages) 
We first define the four-dimensional manifold M with abstract indices {a, /3, ^, v, A, ct}: 

In[2] DefManifold[ M, 4, {a, /3, ^, v, A, cr} ] ; 
and then we define the ambient metric g of negative signature, along with its associated covariant derivative CD: 

In[3] :^ DefMetric[ -1, g[-Q;,-/3], CD, {";","V"}, PriiitAs->"g" ]; 

where M and g respectively correspond to M. and g. Several tensors related to this metric are automatically defined 
at the same time (e.g., all the curvature tensors). Note that in xAct, the covariant indices of a tensor are represented 
with a minus sign (gC-a, -/3] means (jap), while the latter is omitted for eontravariant indices (gCa, /3] means g"^^)- 
Upon defining the perturbations dg of the metric g with the command: 

In[4] ■'= Def MetricPerturbation [ g, dg, e ]; 
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where s is the perturbativc parameter to be used in the expansions, it becomes feasible to evaluate the perturbation 
of any tensor associated with the metric. For instance, the perturbation at first order of the Ricci scalar is simply 
obtained by evaluating: 

In[5] := ExpandPerturbationSPerturbed [ RicciScalarCD [] , 1 ] // ContractMetric// ToCajionical 

Out[5] := i?[V] - e dgi"^i?[V]a0 + e V pVc^dg^"'^ - e VpV^dg^""^ 

The functions ExpandPerturbation and Perturbed is used to evaluate the perturbation of any expression up to any 
order (here at first order), the function ContractMetric removes the background metric tensor through contraction 
on dummy indices (i.e. repeated indices), and the function ToCanonical simplifies the result, gathering together the 
terms which are equal up to symmetries. Further details can be found in Q. 

C. Conformal transformation 

In xAct, the first metric defined on the manifold is the one that is used to raise and lower any tensor indices. For 
our purpose, we have chosen it to be the conformal metric g (cf. input In[3]), which is different from, but conformally 
related to, the background metric of the physical space-time. This choice ensures that the conventional way of moving 
the indices of perturbed fields (equation ([2])) is well satisfied within our algorithm. We however need, now, to relate 
the tensorial quantities one considers in perturbation theory, namely those living on the background manifold of the 
physical space-time, to those we have defined or shall define on the conformal background manifold M. 

We detail this important point in the rest of this section. 

1. Mathematical framework 

Let us denote by g the metric of the physical space-time, and by g its background value. The metrics g and g are 
related by the conformal transformation 

= [g-'r = ^-^{9-'Y\ ig-Tgp.^s^,, (12) 

with a being the scale factor of the background space-timed. Substituting into the first expression of (fT2|) the 
perturbative expansion ([1]) for g and its counterpart for g, we extend the conformal transformation to the background: 

9>.. = a'-9^., Q-'y^a-^r^ (13) 

and to the perturbed level: 

^" V = , with f">^'''' = 5^^" V ■ (14) 

The associated Levi-Civita connections V and V, on the one hand, and V and V, on the other, are related bjH 

V^,uj^ = V^w^ - CP^^ujp , VpUJt, = V^w^ - CPp^ujp , (15) 
for any one- form u. Using equations ([T2|) and (fT3|). we can write the quantities C^^, and C''^^^ as 

CP^^^25\\I,)\na-gp,Vnna, CP ^ 25\\I ,)\na - g^^nna , (16) 



^ We stress again that the first metric defined on M is the one that is used to raise and lower any tensor indices. Defining another metric, 
say /, on the same manifold from the function DefMetric actually creates two objects: (i) the tensor /, with internal notation f ; and 
(ii) the tensor with internal notation Inv[f] . These have the following properties: 

which explains the notation we use in equation l|12|l . 
* By definition the scale factor used in the transformation II12I I is not to be perturbed. Hence we have: g = g, and thus: V = V. 
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with V = gP'^Wa and V = g'"'V a, and where the parentheses indicate symmetrization over the indices enclosed. 
We can then formulate the n**^ order perturbation of a^: 

» I 

With the help of the two previous relations, we are now able to provide the correspondence we seek. 
For instance, the Riemann tensor associated with the metric g is given by (see appendix D of [l6| ) 

= R^,,; + 2 ^'^[^ V,] Vp In a - 2 .gp[^ V,] V" In a 

- 2 (5'"[p V^] In a Vp In a + 2 gp[pSv\ In a V'" In a - 2 5p[p'5''^] V« In a Vc In a , (18) 

where the brackets indicate anti-symmetrization over the indices enclosed. Perturbing this expression and using 

equation ([TT]), we can finally relate A" [i?^^^''] to A" [i?^^^''] and recover the usual quantities studied in perturbation 
theory. 



2.. Implementation in xPand 

The xAct package provides the tools to define a metric conformally related to another, thanks to the option 
Conf ormalTo of the function DefMetric. We have encapsulated this in xPand within the function Def Conf ormal- 
Metric, which furthermore ensures the transitivity of several conformal transformations. 

Let us load the package xPand: 

In[6] := «xAct'xPand' 



Package xAct'xPand' version 0.4.0, {2013,02,08} 

CopyRight (C) 2012-2013, Cyril Pitrou, Xavier Roy and Obinna Umeh under the GPL. 

By evaluating the command 

In[7] := Def Conf ormalMetric [ g, a ] ; 

we define the scalar factor a[] and the metric ga2 conformally related to g through a (ga2 thus corresponds to the 
background metric g of the physical space-time) . To obtain the expression of any tensorial quantities living on the 
manifold described by ga2 in terms of those defined on M. one then simply has to use the xPand function Conformal. 
For instance, we obtain for the Riemann tensor associated with ga2 the expression 

In[8] := Conformal[ g, ga2 ][ RiemajinCD[-Q!,-/3,-/i, i/] ] 

Out[8] := R [VJ - 5p + <5„ gp^ +2 6^ 2 6^ 

-5a hd„ 2g^^ 2 \-^9a^L 2 ^9Pt^ Sap 



^ The background and perturbed connections, on the one hand, and the background and conformal connections, on the other hand, are 
respectively related by 

for any one-form uj. Performing a conformal transformation on the former expression and perturbing the latter respectively yields: 

^p^u = Vpi^^ - CP^^Up = \7pZ^ - A[TP^^] Up - CP^^Up . 
These relations can only be compatible if 

It can be chocked directly from Eqs. (|6]l, I I14I I and 1171 1 that this is indeed the case. This shows the equivalence between the transformations 

V — >■ V — > V and V — > V = V — )■ V. The latter approach actually proves itself to be faster within xPand. It is therefore the one that 
we have coded in the function ToxPand (see section [IV Bl l. 
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which coincides with equation p8|) . 

The conformal transformation may as well be performed on quantities that are not related to a metric. For a 
general tensor, it is defined as 

j,Hl...p,p _ ^q-p+W{T) rptli...flp ^g'j 

where W{T) is the conformal weight of the tensor T. The default value of W is chosen to be zero in order to leave 
the norm invariant under a conformal transformation: 

In[9] := Def Tensor [ W [-a] , M ] ; 

In[10] :== Conformal[ g, ga2 ][ W[-a] ] 

Out[10] aWo, 

In[ll] Conformal [ g, ga2 ] [ WCa] ] 

Outfll] := 

a 

The conformal weight can however be modified for each tensor with the xPand function Conf ormalWeight. This can 
be of use for instance to preserve the geodesic character of light-like vector^: 

In[12] := Def Tensor [ kC-a] , M ] ; 
In[13] Conf ormalWeight [ k ] ~ = -1; 
In[14] := Conf ormalWeight [ k [-a] ] 

OutfU] := 

In[15] := ConformalWeight[ k[a] ] 

Out[15] := -2 

In[16] := ConformalC g, ga2 ][ kCa] CD [-a] Ok [-/3] ] 

Out[16]:^ 

So far, by applying Conformal then ExpandPerturbationSPerturbed on a given expression defined on M, one 
obtains the perturbation of its conformal transformation in terms of the tensors defined on M, the metric g, its 
perturbations dg, the connection V, and the scale factor a. To end up with the usual expressions of perturbation 
theory, one needs to perform a 3 + 1 splitting of the background manifold (section |ll]), then decompose each perturbed 
fields into its spatial and temporal parts and finally parameterize the perturbations of the metric (section IIII[) . 

II. 3+1 SPLITTING OF THE BACKGROUND MANIFOLD 
A. Induced Metric 

The assumption that the background space-time possesses a set of (three-dimensional) homogeneous surfaces pro- 
vides a natural choice for the 3-1-1 slicing. We foliate the background manifold by means of this family, and we denote 
by n the unit time-like vector (n'^n^ — —1) normal to it. The metric of Ai is decomposed as 

g^iu = hf^^i, - n^fiy , with h^uyfi^ = and h^^ph''^ = /i^^, , (20) 

where h represents the induced metric of the spatial hypersufacefl 



^ We will make use of this prescription in a future version of xPand to implement the derivation of the (perturbed) null geodesic equation. 
^ For the sake of clarity, let us note that ^"^h are not the perturbations of h. From the definition of '"'h together with relation I I20II . we 
can actually relate them as: A"[h.] =^"l/i — A"[fi ® n\. 
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The acceleration of the so-called Eulerian observers satisfies in all generality [l2| : 

a^=nP\7pn^^^, (21) 
a 

with a being the lapse function. D stands for the connection of the three-surfaces associated with h {Dph^,^ = 0), 
and it is related to the four-covariant derivative as 

for any spatial tensor fielc0. Since the lapse is homogeneous in the configuration at stake, the acceleration vanishes 
and the observers are in geodesic motion. We can therefore label each hypersurface by their proper time j] and write: 
Tip = — V^?7. In addition, n being hypersurface-forming by construction, its vorticity vanishes; this property yields: 

Cj^, = hP^h%^[pn,] =0 ^ V[,,7-i„] = , (23) 

where the equivalence stems from the null acceleration. For comprehensive reviews on the 3-1-1 formalism, we refer 
the reader to, e.g., [H, 01- 



B. Extrinsic Curvature 



Another tensor we shall make use of is the symmetric extrinsic curvature tensor, which characterizes the way the 
three-surfaces are embedded into the background manifold. It satisfies the relation: 

K^, = hPji%Vpn, , (24) 

where we have chosen a positive sign for the right-hand sid£0. From the decomposition (|20p along with the vanishing 
of the acceleration a and the unitary of n, we can reformulate expression (j24p as 

= Vpfi^ . (25) 

Since the volume expansion of the background space-time is entirely contained in the scale factor a, and owing to the 
conformal transformation (jl3p . the trace of the extrinsic curvature vanishes: K^^ = 0. As a result, we have for general 
Bianchi cosmologies: K^i, = (7^^, with a^^ being the shear of the Eulerian observers; and for FLRW cosmologies: 
Kp, = 0. 



C. Curvature tensors 



The splitting of the four-Riemann tensor can be constructed from its different projections onto the spatial slices 
and the congruence of the observers. It is written: 

Rpypa = ^Rp.vpa + 2/?^[pA'^]^ - 4 (_D[^^^][p) fi^] - 4 (i5[pA'^][p) fi^] + 4 fi[^ K^^ K(_y n„] + 4 K^][pn„] , (26) 

where ^Rpvpa stands for the three-Riemann curvature of the hypersurfaces. The over-dot indicates the covariant 
derivative along the world-lines of the observers (for any tensor field T, we have: T^j...^ = n^'^ pT^j^...p ). The purely 
spatial projection of (|26p only calls upon the first two terms, and it drives the Gauss relation: 

h'Pji\h^ph^, R^^^C - '^A.-P- + Si^MP^-]^' ■ (27) 
The three-space and one-time projection gives, from the next two terms, the Codazzi relation: 

h'^^h^h^^ni i?^,^^ = D^K.p - D.Kpp , (28) 



® We recall that the operator D loses its character of derivative when it is applied to non-spatial tensors. More precisely, we are not 
allowed to use the Leibniz rule anymore, as one can realize upon writing for instance: 

5p(V'f^i...Mj = V£'pT^i...Mp + ^"Vi ■ ■ ■ h"%,T,,...,^D,i, 

for any scalar field '0. 

® This convention does not affect the 3 -|- 1 Einstein equations as written in terms of the kinematical quantities of the Eulerian observers. 
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and the last non-null projection (two-space and two-time) provides, from the last two terms, an evolution equation 
for the extrinsic curvature. These decompositions can be performed in xAct with the function GaussCodazzi. 
For FLRW space-times, the curvature tensors of the hypersurfaces cast the form: 

R^iupij — 2fc ^p[/j ^iy](T : Rp.y — 2A: Jl^i/ , R = 6fc , (29) 

with k being the curvature parameter (equal to zero for flat FLRW cosmologies). The corresponding expressions for 
general Bianchi space-times are more involved, as they require the introduction of the constants of structures. We 
detail their derivation in appendix |B] 



D. Derivatives 

In order to achieve the 3-1-1 splitting, we are left with the decomposition of the covariant derivative V in terms of 
the induced derivative D. For general spatial tensors (namely, for spatial tensors defined within M. or defined within 
M then mapped onto A^), the relation between the two derivatives reads: 

p 

^ pT^i...^p = — ^ip'7)ii...pp + D pT^j^-^ p^^ + ^ ^ K p T'pj...p;_j(jpi+i...pp ■ (30) 

i=l 

Even though our formalism is purely geometrical, we aim at eventually providing, for the perturbations, partial 
differential equations with respect to the proper time rj of the Eulerian observers. When considering the four- 
dimensional basis built to address the Bianchi classification (refer to appendix [B| , the Lie derivative along the 
direction of n precisely comes down to 9,, . It is accordingly more appropriate for our purpose to use the Lie derivative 
rather than the dot derivative. 

The relation between Cn and n'^Vp is writterF°l: 

p 

'^fi^pi...pp = 'j)ii...pp + p; 7pj...p._j(jpi+i...pp I (31) 

1=1 

and provides us with the following reformulation of (j30[) : 

p 

^ pTp-^^...pp ~ —fipCfiTp-^ p^ + DpTp-^ p^ + 2 ^ ^ fl(^. K p) Tpj...p._jcrpi+i...pp • (32) 

i=l 

Finally, a last expression we shall need is the commutation rule between the derivatives Cn and D. For general spatial 
tensors, it is given by 

Cn {DpTp„„p^) = Dp {CnTp,...p^) + ( /i^'^^c^^pm. " ^p^^p," - ^m.^^p' ) T^mi-m.-i-m.+i.-.m. , (33) 

i=l 

where we have made use of relations ((22l) and (|28p for its derivation. 



E. Implementation in xPand 

The 3+1 splitting of the background manifold is performed by the xPand function SetSlicing. It can be applied 
to the following spatially homogeneous cosmologies: "Minkowski", "FLFlat", "FLCurved", "Biancliil", "BianchiA", 
"BicinchiB" and "Anisotropic". 



Note that for a spatial tensor T, the quantity CnTp^,,,^^ is also spatial. 
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1. Construction of the spatial hypersurfaces 

From the ambient metric g, SetSlicing first defines the unit normal vector n, the induced metric h of the hy- 
persurfaces, the associated covariant derivative cd, the associated scale factor a[h] , and the conformal metric gah2. 
It then specifies the expressions of the intrinsic and extrinsic curvature tensors according to the type of cosmologies 
chosen by the user (hence, the geometry of the model is fully described). 

Let us illustrate with the cosmologies "FLCurved" and "BianchiA" (for the latter most sophisticated case, refer to 
appendix IB]) . With the evaluation 

In[17] := SetSlicingC g, n, h, cd, {"|","D"}, "FLCurved" ]; 

the extrinsic curvature K [h] [-a , -/?] is set to zero, and the expressions of the three-curvature tensors arc implemented 
following relations ([29l) . For instance, we have for the three- Ricci tensor: 

In[18] := Riccicd[-Q;,-;3] 

Out[18] := 2hai)k 

For "BianchiA" cosmologies, 

In[19j := SetSlicingC g, nA, M, cdA, {"T'/'D"}, "BianchiA" ]; 

SetSlicing defines the spatial constants of structure C'^ij (see appendix IB|) and then constructs a function allowing 
one to express the intrinsic curvature tensors in terms of them. For the threc-Ricci tensor, we then have: 

In[20j := RiccicdA[-b,-c] // ToConstantsOf Structure [hA] 



Out[20j 




which is equivalent to equation (|B2ip F^ . The constants of structure can be further expanded, following the usual 
Schiicking, Kundt and Behr (SKB) decomposition [l7| . by means of the xPand function ToBianchiType [M] . 

We summarize in table|T]the different evaluations performed by SetSlicing with respect to the type of cosmologies. 



Space-time 


Extrinsic curvature 


Three- curvature tensors 


Constants of structure 


Minkowski 


Null 


Null 


Null 


FLFlat 


Null 


Null 


Null 


FLCurved 


Null 


Eqs. 


C- ij — ^y/k h ' €mij 


Bianchil 


K^, 


Null 


Null 


BianchiA 




Eas. (|B20|). (|B21|). (|B22|) 




BianchiB 


K^, 


Eas. (|B20(). (|B21|). (|B22|) 




Anisotropic 




Eas. (|B20|), (|B21|). (|B22|) 


^ ij 



TABLE L Evaluations performed by SetSlicing for the extrinsic curvature, the three-curvature tensors and the constants of 
structure, according to the type of homogeneous cosmologies. In terms of g (or, equivalently, g), the dynamics of "Minkowski" 
and "FLFlat" space-times are identical. The difference lies in the conformal transformation (|13|) : for the former models, the 
scale factor is set to one. For "BianchiA", "BianchiB" and "Anisotropic" space-times, the three-curvature tensors can be 
formulated in terms of the constants of structure by means of the function ToConstantsOf Structure [] and can be further 
expanded with the help of ToBianchiType [] . The "Anisotropic" space-time is used for hypersurfaces whose dimension differs 
from three (for these models the SKB decomposition does not apply). 



Equation IIB21I I is recovered by making use of the Jacobi identity IIB5I I on Out[20j. This latter relation is not implemented, as it implies 
a symmetry among several terms that xAct does not yet handle. However, once the constants of structure are expanded according to 
the parameterization IIB8I I. the Jacobi identity reduces to the constraint I IBQI l. which is automatically applied in xPand. 
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2. Derivatives 

For anisotropic cosmologies, SetSlicing creates a set of automatic rules to handle the action of the derivatives 
on the three-curvature tensors. The aim is to obtain a formulation in terms of the extrinsic curvature tensor and/or 
the constants of structure, which fully identify the geometry of the model. Note that such rules are not necessary 
for isotropic models: the expressions of the intrinsic curvature tensors are simple enough to let xAct perform the 
evaluation. 

The action of the derivative V is implemented following equation (|32p . that of D following equation (jB14[) 
together with (jB16p . and the application of Cn on a three-curvature tensor automatically calls the function 
ToConstcintsOf Structure [] so as to use the property (|B19p . Regarding the three-covariant derivative, we have 
for instance: 

In[21] := cdA[-Q;]@RiccicdA[-^,-i/] 

Out[21] \ Cj R{b\,p - \ R{DU + \ Cj R{b\,p -f \ C„/ R\I)U " \ R[bU 

Lastly, SetSlicing constructs an automatic rule to perform the commutation p3p for any expression. This serves 
to make sure that the Lie derivative will first act on a given tensor, so as to recover the usual formulation of the 
perturbed equationj^. For "FLCurved" cosmologies, we have: 

In[22] := Def Tensor [ VC-a], M, OrthogonalTo->{n[a] } , ProjectedWith->{h[a,-/3] } ]; 
In[23] :^ LieD [n [/x] ] @cd [-i^] OV [-a] 
Out[23] := Cn Va 
and for "BianchiA": 

lnl2^] := Def Tensor [TA[-a!],M,OrthogonalTo->{nA[a] },ProjectedWith->{hA[a,-/3] },PrintAs->"T"] ; 
Inl25] := LieD [nA ] OcdA [-i/] OTA [-a] 

where to obtain Out[25], expression ()B14p has been automatically applied to the extrinsic curvature tensor. 

III. PERTURBED FIELDS 

We now present the splitting of the perturbed fields into their spatial and temporal parts, the scalar- vector-tensor 
(SVT) decomposition of their projected components, and finally the parameterization of the metric perturbations. 

A. Projected components and scalar-vector-tensor decomposition 

Any perturbed quantity lying within and mapped onto the background manifold can be decomposed by means 
of the normal vector n and the induced metric h. For a rank-two covariant tensor T, we have for instance: 

T^. = n^n,{nPn-'Tp„) + 2ni^^{nPh\)Tp„) + {hPji%Tpa) . (34) 

The terms inside parentheses hence constructed respectively define a scalar, a spatial vector and a spatial tensor. As 
is customary in perturbation theory, we shall then employ the SVT decomposition to further expand the two last 



In line of this comment, let us note that SetSlicing also creates internal rules for the commutation of several D. These will enforce 
the appearance of Laplacians, and move a covariant derivative closer to a tensor when a divergence is present. 
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projected quantities [l8|. Let us briefly recall the form of this decomposition for later needs: a spatial vector [7^ is 
split into a scalar part S and a vector part as 

Uf, = D^S + V^, with L)^y^ = 0, (35) 

and a symmetric spatial tensor i?^^ is decomposed into two scalar parts 5*1 and 15*2, a vector part VJ^ and a tensor 
part Tpi^ as 

H^. = D^D,Si+S2h^.. + D^^V,)+T^., with l^l^'^^J^^ ^^^^ (36) 



B. Perturbations of the metric 



The SVT decomposition of the metric perturbations yields the general expressions: 

n^n-f-V = -2^'"'/', (37) 

nPh%'"X, = -D, - f">B, , (38) 

h'^h", '"V = 2 (^^^, t">£; + D^^ + f">S^, - f"V ^M-) • (39) 

Four of the ten degrees of freedom of are carried by the scalars *"^(/), and ^"*_B; four are encoded in the 

vectors ^"^E^ and ^"^5,^; and two are contained in the tensor ^"^''E^i,. Note that some of these fields are required to 
vanish for specific gauge choices (see, e.g., [l^ for a comprehensive review). 



C. Implementation in xPand 

1. Construction of SVT quantities 

The construction of spatial tensors satisfying the SVT properties is performed by the xPand function Def Pro j ected- 
Tensor. By default, these tensors are defined on both the background and perturbed manifolds. For instance, we 
have: 

In[26] :^ DefProjectedTensor[ U[-a,-/3] , h ]; 
In[27] n[-a] U[a, /3] 

Out[27] :^ 
In[28] := U[a, -a] 
Out[28] 

In[29] := Perturbation [ U [-a, -^] , 1 ] 

Out[29] A[C/a/3] 

To relax one or several default properties, the user has to modify the optional arguments TensorProperties, set by 
default to {"SymmetricTensor" , "Traceless" , "Transverse"}, and SpaceTimesDf Def inition, evaluated by default 
as {"Background" , "Perturbed"}: 

In[30] Undef Tensor [ U ] ; 

InfSl] DefProjectedTensor [U[-a,-/3], h, TensorProperties->{"SymmetricTensor" , "Transverse"}, 
SpaceTimesDf Def inition->{ "Background"}] ; 

In[32j := U[a,-a] 

Out[32] :^ U\ 

In[33] := Perturbation [U [-a, , 1] 

Out[33] := 
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2. Comments on the label-indices 

In the internal notation, xPert attaches a label-index to the metric perturbations in order to denote their order. 
The first order perturbation dg^^ (cf. output Out[5]) is hence stored as dg[LI [1] ,-a,-/3] . In xPand, we employ the 
same notation for spatial tensors defined with Def ProjectedTensor, and we moreover attach a second label-index to 
indicate the number of Lie derivatives along n acting on them. Thus, we have for instance: 

In[34] := DefProjectedTensor [V[-q:] , h] 
In[35] := LieD [n[/z] ] (§V [LI [1] , LI [0] , -a] 

Out[35] :^ 

In[36] :== V[LI[1], LI[1], -a] 

Out[36] (I'K^' 

DefProjectedTensor constructs a set of rules to automatically allocate label-indices to a tensor written without. 
Hence, V[-q;] is converted to V[LI[0] ,LI[0] ,-q;], while V[LI[p] ,-a] is converted into V[LI[p] ,LI[0] ,-«], for 
any perturbation order p. 

Note that xAct interprets a tensor with label-indices as a tensor by itself, whatever the meaning of the label-indices. 
Its indices arc therefore raised and lowered by the ambient metric in the usual way: gCa,/?] V[LI[1] ,LI[1] ,-a] 
yields V[LI[1] ,LI[1] ,/3]. While this is obviously mathematically correct for isotropic manifoldj^. it is no longer 
true for anisotropic ones. For these latter models, the second label-index can be interpreted as Lie derivatives only 
when the tensor is with indices down. To avoid confusion, we modify the output when the tensor indices arc up. For 
instance, we have for "BianchiA" cosmologies: 

In[37] := DefProjectedTensor[VA[-a] , M] 

In[38] := VA[LI[1], LI[1], -a] 

Out[38j ~ (I'K,' 

In[39] := VA[LI[1], LI [1] , a] 

Out[39] 

3. Construction of the perturbations of the metric 

The perturbed fields introduced in subsection IIII Bl are constructed with the xPand function Def MetricFields. 
The evaluation of the command 

In[40] := Def MetricFields [ g, dg, h ]; 

calls the function DefProjectedTensor in order to define the projected components of the metric perturbations and 
allocate them all the SVT properties. The set of rules (|37)). (|38|) and (pQ)) arc then automatically defined using the 
xPand function SplitMetric: 

In[41] := GaugeRules = SplitMetric [ g, dg, h, "AnyGauge" ] 
In[42] := dg[LI[l], -a, /. GaugeRules 

Out[42j :^ 2 (I'i^a^ + ^^'>Bf) ~ (^'^^ fi/j - 2 n„ - 2 - np DJ^^B + (■^''Ef, - n„ Dfj (^'S 

The particular gauges we have implemented are: "ComovingGauge", "FlatGauge", "IsoDensityGauge", "Newton- 
Gauge" and "SynchronousGauge" . Should the user wish to consider other gauges, own rules can be created as follows: 
MyGauge = {dg[LI[ordJ, ^_ :> ...}. 



In such cases, g[a,/3] V[LI[1] ,LI[1] ,-a] corresponds to g°'^ CnVa which is equal to CaV^ , as the extrinsic curvature vanishes. This 
last expression indeed corresponds to V[LI[1] ,LI[1] 
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4- Splitting of the background covariant derivative 

We now have nearly all the necessary tools to obtain the final expression of any perturbed field. Let us quickly 
review the previous steps and then introduce the last function we need. 

To derive the perturbation of, e.g., the four-dimensional Ricci scalar, we first perforin a conformal transformation 
from g to gah2 (and express it in terms of the metric g), then we perturb the resulting expression, and finally we 
substitute the metric perturbations by their SVT components using the set of rules GaugeRules: 

In[43] := Conf ormalCg, gah2] [ RicciScalarCD [] ] 
In[44] := MyR = ExpandPerturbationOPerturbed [ 7., 1 ] 

In[45] := (MyR /.GaugeRules) // ProjectorToMetric // GradNormalToExtrinsicK 
// ContractMetric // ToCanonical 

The final result has been expanded using the xAct functions ProjectorToMetric, which replaces /i^^ by g^^ — fif^n^, 
and GradNormalToExtrinsicK, which replaces V^n^ by K^j^^ and it has been simplified with ContractMetric and 
ToCanonical. 

Since the relation p2p is not automatically evaluated so far for quantities other than the three-curvature tensors, the 
result of In[45] still involves the covariant derivative V. In order to split the latter in terms of the induced derivative 
and Lie derivative along n, one finally needs to use the xPand function SplitPerturbations, which applies as well 
the Gauss-Codazzi decompositions. In the above example, we obtain: 

In[46] :^ SplitPerturbations [7, ah [] '2 , h] 

Out[46] := 6H^ + 6H' + 6K: + e{-12-H^ '^'0 - 12 W '^'0 - 6 H + 12 IC-ISH <^y - 6 ''V" 

+6 H DaD°' 'I'B - 2 DcD°' '^'B' + 6^. 0^0°"-^^' + 2 <i'i;" - 2 <i'0 4- 4 D^D^'^V) 

This way of proceeding is, however, rather inefficient for general gauges at higher order. Instead, the set of rules 
GaugeRules can be used in an optimized manner by the function SplitPerturbations itself: 

In[47] := SplitPerturbations [ah [] 2 MyR, GaugeRules, h] 

In such a way, the rule following relation ([32]) is evaluated on the projected components of the metric before specifying 
their SVT decomposition, which takes much less time at higher orders. 

IV. FEATURES OF THE ALGORITHM 
A. Summary 

Let us review the main steps that need to be followed in order to derive the perturbation of any fields. 

i. We define the background manifold M (with DefManif old) and the ambient metric g (with DefMetric). The 
splitting of M is realized with SetSlicing, according to the type of cosmology chosen by the user. 

ii. We apply to the quantity to be perturbed a conformal transformation from g to g and express the result with 
respect to the quantities defined on M. This is done with the function Conformal. 

iii. We use the xPert tools Perturbed and ExpandPerturbation to perturb the previous expression at any order, 
in terms of the metric perturbations dg and other tensors. 

iv. We use the functions Def MetricFields and SplitMetric to define and construct the SVT parameterization of 
the metric perturbations (thanks to the function Def TensorProjected), according to the gauge chosen by the 
user. The perturbations of the fluid quantities are defined and decomposed with the functions Def MatterFields 
and SplitMatter (see further below). 

V. We finally use the function SplitPerturbations to decompose the covariant derivative V in terms of the 
induced derivative and the Lie derivative along n. The Gauss-Codazzi relations are also applied, and the 
constructed rules are used to transform the resulting expression. 

These four last steps have been coded all at once in a single function called ToxPand. We present it further below, 
after introducing a minimal example. 
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B. A minimal example 

We now propose a brief and self-contained example to illustrate our package. The parameterization of the metric 
perturbations is here constructed by hand, and we only consider the Bardeen potentials. 

<<xAct/xPand.in; 

DefManifold[ M, 4, {a, (3, fi, v} 'i; 

DefMetric[ -1, g[-a,-^], CD, {";", "V"} ]; 

DefMetricPerturbation[ g, dg, e ] ; 

SetSlicingC g, n, h, cd, {"I", "D"}, "FLCurved" ]; 

order = 1 ; 

Def ProjectedTensor [ (/)[], h ]; 
Def ProjectedTensor [ ^jLI , h ]; 

MyRicciScalar = ExpandPerturbationOPerturbed [ Conformal[g, gah2] [RicciScalarCD [] ] , order ]; 
MyGauge = dg[LI[ordJ, a_, :> - 2 n[a]n[/3] (^[LI[ord]] - 2 h[a, /3] VELlCord]]; 

SplitPerturbations [ ali[]~2 MyRicciScalar, MyGauge, h ] 

The output generated by the last line is (compare with Out[46j) 

Out[47] 6 + 6 ^' + 6 /C + e (-12 d,^ _ ^2 W <^'(/) - 6 H <^y + 12 '^V/C - 18 H <^V' - 6 
-2 DaD"' '1)0 + 4 (^V) 



C. Secondary functions 

Even though the user is free to parameterize tensor perturbations by creating its own rules with projected tensors 
defined from Def ProjectedTensor, we have seen in paragraph IIII C"3l that the functions Def MetricFields and 
SplitMetric can take care of this procedure for the metric perturbations. Similarly, we have implemented in xPand 
the functions Def MatterFields and SplitMatter for the definition and parameterization of the fiuid perturbations, 
and more precisely for those of the energy density, pressure and fiuid four- velocity. 

We also wish to mention that we have extended in xPand the xPert function GaugeChange which performs gauge 
transformations at any order for a given expression. Our extension SplitGaugeChange executes a 3 + 1 splitting of 
these transformation rules. For a glance at the gauge transformations for metric perturbations and the construction 
of gauge-invariant variables, see [lH, [l^ [l^, [2l|. For examples of use of the secondary functions SplitMatter and 
SplitGaugeChange, we invite the reader to go through the example notebooks which are distributed along with the 
package xPand. 

Furthermore, for the cases where one wants to consider the predefined gauges, it is enough, simple and straightfor- 
ward to use the xPand function ToxPand to obtain the perturbation of any expression. For instance, the following 
five lines suffice to derive that of the four-Ricci scalar in any gauge and up to order two: 

<<xAct /xPand . m ; 

DefManifold[ M, 4, {a, (3, fi, v} 1 ; 
DefMetric[ -1, g[-a, -/3] , CD, {";", "V"} ]; 
SetSlicing[ g, n, h, cd, {"I", "D"}, "FLCurved" ]; 
ToxPand[ RicciScalarCD [] , dg, u, du, h, "AnyGauge", 2 ] 

where u and du are the fiuid four-velocity and its perturbation, respectively. The function ToxPand combines several 
functions so that the user can easily obtain the desired perturbations in a given gauge without having to deal with 
any detail of the algorithm. 

Finally, the xPand function ExtractComponents allows to extract the projected components of any tensor. For a 
given background slicing, the user only needs to specify the type of projection (see appendix [X] for an example). Note 
also that for a rank-one (resp. rank-two) tensor, the xPand function VisualizeTensor allows to display all projected 
components in a vector (resp. matrix) form. Again, we refer to the example notebooks which are distributed along 
with xPand, for more details about the syntax of these functions. 
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D. Recovering standard results 

We have checked that with our implementation we recover the standard resuhs of cosmological perturbation theory. 
More precisely, our algorithm is in accordance with 

• all first-order results of Einstein equation and stress-energy tensor conservation equation, for flat and curved 
FLRW space-times in any gauge (see for instance Ref. p^): 

• all first-order results for Bianchi type I background, in the gauge chosen in Refs. [22l. [23j: 

• all second-order results of Einstein equation and stress-energy tensor conservation equation, for flat and curved 
FLRW space-times in the Newtonian gauge (see for instance Ref. [l^ for the complete set of equations). 

Our package now enables to extend these results to higher order and in any gauge. It also allows for the study of any 
type of perturbed Bianchi cosmologies. 

It is worth mentioning that in order to obtain useful standard differential equations with respect to the conformal 
time 77, it is necessary to perform a mode expansion on the hypersurfaces, that is one needs to find the eigenmodes 
of the spatial Laplacian D^D^. This is simple for flat FLRW cosmologies, where one just has to employ a Fourier 
transformation, and it is also well-known for curved FLRW models, where hyper-spherical Bessel functions need to 
be used [131 ■ However, the procedure is still unknown for general Bianchi cosmologies. Apart from the special case 
of Bianchi type I, where the modes can also be found from a Fourier transformation (and thus lead to a simple set of 
equations |25|). there is no general technique to obtain the eigenmodes of the Laplacian for all other types. Only in 
special cases (see for instance Ref. [1^) this has been done explicitly. 

E. Timings 

In practice, the timing for the computation grows like (slightly faster than) a power law of the perturbation order, 
whatever the gauge (see Fig. [1]). It takes xPand less than two minutes (see Fig. [2]) to decompose completely the 
perturbation of a rank-two curvature tensor, such as the Ricci or Einstein tensor, up to second order in any gauge. 
The decomposition for the perturbations of the Riemann or Weyl tensor, up to second order and in any gauge, takes 
a little more than two and thirteen minutes, respectively. 

CONCLUSION 

xPand is the first comprehensive package that allows to perform algebraic calculations in cosmological perturbation 
theory for homogeneous background space-times, up to any order and in any gauge. It is worth stressing again the 
several features and advantages of this package. 

• The expression of any perturbed field can be obtained for all Bianchi cosmologies, up to any order of interest 
and in any gauge, in a very simple and straightforward way. 

• The package avoids the complexities of a component-by-component computation of the perturbed fields thanks 
to the use of the 3-1-1 formalism. 

• The package is relatively fast: at first order, all the perturbed equations can be derived in approximately two 
seconds, and as the order of the perturbation increases, the timing grows roughly as a power law. 

• The computations can be applied to space-times of arbitrary dimension, and for any extension of general 
relativity (the current exception being gravitational theory with torsion) . 

• Taking advantage of the package xTensor, xPand handles tensor indices just the same way a user would do when 
calculating by hand. In particular, it does not break the summation over repeated indices (as e.g. in Out[4-6]). 

• It totally eliminates the laborious summation over repeated indices, an element that has inhibited the use of 
other packages developed to solve similar models. 

• The user who has little knowledge of the Mathematica or xTensor syntax can obtain from the package almost 
the same utility as an expert would do. 

Finally, we plan to extend the scope of xPand beyond the derivation of perturbed equations in general relativity, 
so as to provide the entire Einstcin-Boltzmann system for radiation transfer, and the Einstein-Jacobi map system 
needed for understanding the effect of weak gravitational Icnsing. 
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Perturbation of Ricci scalar in synchronous gauge Perturbation of Ricci scalar in flat gauge 
10' [ 1 10' [ 




1 2 3 4 1 2 3 4 

Order Order 

FIG. 1. Timings for the perturbation of the four-dimensional Ricci scalar in different gauges. From top to bottom and left to 
right: synchronous gauge, spatially flat gauge, comoving gauge, and Newtonian gauge. On each plot, the curves from bottom to 
top refer to: (i) formal perturbations with xPert using a conformal transformation (red line); (ii) perturbations for a Minkowski 
background (yellow); (iii) perturbations for a curved FLRW background (green); (iv) perturbations for a Bianchi I background 
(blue); and (v) perturbations for a general Bianchi background ("BianchiB") (purple). All timings were performed on a single 
4 GHz core, with a 8 GB RAM. 



Curvature fields in Newtonian gauge Curvature fields in Synchronous gauge 




1 2 3 4 1 2 3 4 

Order Order 



FIG. 2. Timings for the perturbations of the four-curvature tensors, for a curved FLRW space-time, in different gauges. Left: 
Newtonian gauge; right: synchronous gauge. On each plot, the curves from bottom to top refer to; (i) Ricci scalar (red line); 
(ii) Ricci tensor (yellow); (iii) Einstein tensor (green); (iv) Riemann tensor (blue); and (v) Weyl tensor (purple). All timings 
were performed on a single 3.40 GHz core, with a 8 GB RAM. 
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Appendix A: What xPand can do: a non-trivial example 



We here described by way of an example how xPand can be used to derive relativistic perturbed equations for flat 
FLRW backgrounds, up to first order and in the Newtonian gauge. An extension to any other background, order or 
gauge simply requires a change of the related arguments in the following example. 

As usual, we define the geometry of the model by evaluating the following commands: 

In[l] := «xAct/xPand.m; 

DefManifoldC M, 4, {a, (3, fi, v} 1 ; 

In[3] := Def Metric [ -1, gC-^i, -v^ , CD, {";", "CD"} ]; 
Def MetricPerturbation [ g, dg, e ]; 

In[5] := SetSlicing[ g, n, h, cd, {"1", "D"}, "FLFlat" ]; 
order = 1 ; 

For simplicity, we customize the use of the command ToxPand as 
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In[7] := MyToxPand [ expr_, gauge_, order_] : = ToxPand [ expr , dg, u, du, h, gauge, order] ; 

Wc define the constants k, — 8ttG and A with 

In[8] := Def ConstantSymbol [ k ]; 

DefConstantSymbolC A ]; 

and then the energy-momentum tensor of tlie fluid: 

In[10] := DefTensorC T[-^, -ly] , M ] 

We introduce a boolean variable to switch on or off the pressure of the fluid filling the space-time. This allows to 
consider either dust- or radiation-dominated era in a simple way. 

Infll] := $Dust = True; 

In[12] :== IndexSetC T , (pu [] + If [$Dust , , Pu [] ] ) u [a] u [/?] + If [$Dust , , Pu [] ] g[a,/3]]; 

The fields pu (energy density), Pu (pressure) and u (fluid four- velocity) will be automatically created later, when 
calling the function (My) ToxPand. The user is encouraged to set the boolean variable to False in order to include 
the effects of pressure. 

We now derive the perturbations of Einstein equations for the model at stake. We define the Einstein equations as 

In[13] := MyGR[/^_, vj\ := EinsteinCD [/^ , i^] + gLfJ.,!^! A/k - k T[p,i^]; 

and we evaluate their perturbations simply using 

In[14] := MyGRresult = MyToxPand[ MyGR[p, -i^] , "NewtonGauge" , order ]; 

Since the computation of this expression is the time-consuming part, we have stored its value in the variable 
MyGRresult. 

The background components of the resulting expression are extracted using the xPand function ExtractOrder [expr , 
0] (where indicates the background value), and the projected components with the xPand function Extract- 
Components; 

In[15] :^ ExtractComponents [ExtractOrder [ ah[]'2 MyGRresult, ] , h, {"Time" , "Time"}] 

^ r , Aa^ 9 , 
Out[15] := + Ka^p 

K 

In[16j := ExtractComponents [ExtractOrder [ah [] ~2 MyGRresult, 0], h, {"Space" , "Space"}] 

Aa^/jM 

Out[16] := - h^^ -2 h^"^ H' 

K 

Similarly, the first-order equations are obtained by evaluating: 

In[17] := ExtractComponents [ExtractOrder [ ah[]"2 MyGRresult, 1 ] , h, {"Time" , "Time"}] 

Out[17] :^ K (^'p + 6 '^'^ + 6 H '^V' - 20^0° '^'^ 

/n/i^y ;= ExtractComponents [ExtractOrder [ ah[]'2 MyGRresult, 1 ] , h, {"Time" , "Space"}] 

Out[18] a2 (i)B^ p-na^ '^'K p + 5 ^'^B^ - k a^p - 2 HD^ '^'^ -2D^ V' 

/n/iP/ ExtractComponents [ExtractOrder [ ah[]'2 MyGRresult, 1 ] , h, {"Space" , "Time"}] 

Out[19] 2 _ 2 d'S^T^' + Ka^ '"V^^ p + D^ -\D„D'^ '^'S^ + Ka^ pD^" ^ W 

/n/^O/ ExtractComponents [ExtractOrder [ ah[]'2 MyGRresult, 1 ] , h, {"Space" , "Space"}] 

Out[20] ^^^e" + -H + 2 /i^^ + 4 /i^^ H' '^y + 2h>'^U + 4:hi'^n (1 V' + 2 /i^^ V" 
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Appendix B: Bianchi cosmologies 
1. Constructing a four-dimensional basis 

We here review the general properties of Bianchi space-times. The reader can find a detailed presentation of their 



classification in [17|, and summaries in, e.g., j27l l28l|. 

Bianchi space-times posses by definition a set of (three-dimensional) homogeneous hypersurfaces. One can therefore 
introduce three linearly independent spatial Killing vector fields (KVF) ^j, with i S {1,2,3}, satisfying: 

= ^ V(^6.)=0, ri^ef^O. (Bl) 

From these properties together with equation ([23]), we obtain: 

Cn^^^[n,^^]=0. (B2) 

The nature of the Bianchi model is determined by the spatial structure coefficients C'^ij, defined from the commutators 
of the KVF: 

= -C^,a, with C\^=-C%. (B3) 

The Jacobi identity 

[^^, fe, + fe, fc, e.]] + = , (B4) 

contrains these constants to verify: 

C™[yC''j.]„j = => C™'ijC\^ = , (B5) 

where we have used the fact that the C'^ij arc constant on spatial slices. 

We now construct, on a given hypcrsurface, a vector basis {ei} and its dual {e*} invariant under the action of the 
KVF, namely satisfying: 

e, = e, ] = , = . (B6) 

From these properties along with relation (|B3| . we infer: 

[e„ e,] = C%ek , 2e,''e/V[^e^] = -C^ ■ . (B7) 
The constants of structure can be further developed in terms of a symmetric 'tensor' N^^ and a 'vector' A' as 



C% = e.,™^™'^ + 2A[,(5^-,i . (B8) 



where eijm denotes the totally anti-symmetric Levi-Civita symbol. Note that the Jacobi identity (jBSp translates in 
that case to the simple relation: 

N,jA^ = . (B9) 

We then extend the bases {e.;} and {e'} to the whole space-time by Lie dragging them with n, which implies the 
properties: 

Cne^ = [n,e,]=0, Cne' = . (BIO) 

With the above procedure, we are able to construct a four-dimensional basis {sa} = {n, e,} along with its dual 
{e"} = {n, e*} (where n is the dual form of n and a e {0,1,2,3}), that are invariant under the action of the 
KVF. The commutation relations of these new bases simply follow from expressions (|B7|) and (|B10|) : the structure 
coefficients C^^j, vanish when any of the indices is zero and take the values C'^^j otherwise. This method to build a 
four-dimensional basis out of a three-dimensional one defined on a given spatial hypcrsurface is the simplest ontP^. 



Note that in this framework, only fi is a unit vector. An alternative approach consists in building a basis of vectors that all are unitary, 
by renormalizing the e^. However, by doing so, the spatial vectors hence constructed do not commute with fi anymore, and their 
associated structure coefficients become time-dependent. We shall not consider such possibility in the present paper, but details can be 
found in, e.g., [13, mUl. 
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2. Expression of the geometrical tensors 

The components of the induced metric are written with respect to the Bianchi basis {e^} and its dual as 

h^^ = hij e\e\ . (BH) 

From relations (|B1[) and (jB2[) . we deduce: 

£^,V = 0, (B12) 

and using equation (|B11|) together with (jB6|) . we obtain: ek{hij) = 0. The components hij are thus only time- 
dependent, namely hij ~ hij{r]). Any tensor field T lying on the background space-time possesses the same symme- 
tries, and so we can expand it with respect to the Bianchi bases as 

T^,...M. =^.i.....('7)e*Vi---e'V.- (B13) 
Denoting by "^r^fc the coefficients of the connection D in the Bianchi bases, we then deduce that 

p 

i=i 

This relation is used in xPand to compute the induced covariant derivative of any background tensor, such as the 
extrinsic curvature. The connection coefficients associated with the background connection V are given in the bases 
{ea} and {e"} by 

f a&c = ^ ( - ea {(jbc) + et (gca) + Be {{jab) + Cabc ^ Cbca + Ccab) ■ (B15) 

Given that the components hij only depend on 77, we deduce that the spatial-connection coefficients '^Tijk are expressed 
only in terms of the constants of structure. We have indeed: 

^^ijk — ^ijk — l^if^ijk ~ C jki Ckij) , Cijk — ^ijk ~ ^ikj ■ (B16) 

Note that the tensor indices in the bases {ej} and {e'} are lowered with hij and raised with its inverse h^^ , so that 
for instance Ckij = hkmC™ij. Note also that the constants of structure of the basis {ci}, that we here note C[e]'^j^ 
for the sake of clarity, are the components, in this specific basis, of a tensor. The associated components in a general 
basis can be recovered frorrP^ 

C[ef^^ ^ C[e]\^ ek"e\e\ . (B17) 

Similarly, the components of Nij and Ai in a general basis are found from 

7V[e]^, = iV[e]y e'^e^ , A[e]^ = A[e], e\ . (B18) 

Relation (|B14[) is also used to compute the induced covariant derivative of these three tensors, since they also live 
on the background space-time. Note, finally, that from equation (IBIOP and the fact that the C'^j, arc constant, we 
obtain the useful relation: 

Cn{C''^,)=Q. (B19) 

The three- Riemann tensor of the hypersurfaces can be expressed only in terms of the constants of structure. In the 
bases {e^} and {e'}, its components are given by 

- - + IcJfiP", + Cj^C,'" + C^^C^" + -t- ic/,C/^ + ^C'.,C'/ + C^pC/^ , (B20) 



Note that we are not considering the constants of structure of another basis with this construction. We rather build a tensor C[e] = 
Cfe]*. ej. (gi e* (gie-'. It is obvious that the components of this tensor in the basis {e^} and {e'} are the C[e]*'^^, but its components 
can also be taken in a general basis even though we only refer to the commutation structure of the e^. Since xTensor manipulates only 
abstract indices (that is indices in a general basis) , this covariant point of view is necessary to implement the structure of the Bianchi 
space-times in our package. 
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where a double anti-symmctrization [ij] and [kl] is implied on the indices in the right-hand side. The three-Ricci 
tensor and three-Ricci scalar can then be deduced, and we obtain: 

^Rij ~ — -CkiiC^j^ — -CkiiC^j^ + -^Ci'^Cjki + C(ij)^C^pk , (B21) 
^R = - Ja.feC*^*'- - ^a^kC^^'' + C'hCPpj . (B22) 

Note that owing to the Jacobi identity (jBSp . the three-curvature tensors can take several equivalent forms. Finally, 
due to the homogeneity of the hypersurfaces, any induced derivative acting on the three-curvature tensors can be 
computed using Eq. (|BT4)) . All the rules of this appendix, namely equations ((BT4)) . ((BT6| . ((B20| . ((B2T|) and ((B22|) . 
are automatically created in xPand when calling the function SetSlicing, in case the specified space type is a Bianchi 
space-time. The boolean variable $OpenCoiistaiitsOf Structure controls whether or not the constants of structure 
should be opened in the final expressions using the parameterization jBSj. 



Appendix C: Principal commands of xPand 



Conf ormal 


Conf ormalWe ight 


DefMatterFields 


DefMetricFields 


Def ProjectedTensor 


ExtractComponents 


SetSlicing 


SplitGaugeChange 


SplitMatter 


SplitMetric 


SplitPerturbations 


ToBianchiType 


ToConst ant sOf Structure 


ToxPand 


VisualizeTensor 



TABLE II. Information about these commands can be obtained within Mathematica by evaluating ?NameOf Command. 



